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Abstract 

Oh 1 In minimal supersymmetric models the Z-pcnguin usually provides sub-dominant contri- 

butions to charged lepton flavour violating observables. In this study, wc consider the su- 
persymmetric inverse seesaw in which the non-minimal particle content allows for dominant 
contributions of the Z-penguin to several lepton flavour violating observables. In particular, 
and due to the low-scale (TeV) seesaw, the penguin contribution to, for instance, Br(/i — > 3e) 
and fi — e conversion in nuclei, allows to render some of these observables within future sen- 
^\ • sitivity reach. Moreover, we show that in this framework, the Z-penguin exhibits the same 

non-decoupling behaviour which had previously been identified in flavour violating Higgs de- 
cays in the Minimal Supersymmetric Standard Model. 
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1 Introduction 



In recent years, lepton physics has experienced an unprecedented experimental development. A 
non- vanishing - even unexpectedly large - value of the Chooz angle (#13 ~ 9°) has recently been 
measured by several (independent) collaborations p3-[5]. Such a value of #13 opens the door to very 
appealing phenomenological possibilities, among which CP violation in the leptonic sector stands 
as the best example. In parallel to these achievements, the other neutrino oscillation parameters 
(solar and atmospheric) are being determined with very good precision [6j,[T] . In the near future, 
one does expect to identify the fundamental ingredients of flavour violation in the neutral lepton 
sector, but on the other hand, flavour violation in the charged lepton sector still remains to be 
observed. The present and future generations of high-intensity facilities dedicated to discovering 
flavour violation in charged lepton processes render feasible the observation of such an event in 
the near future. 

Despite the fact that minimal extensions of the Standard Model (SM) can easily accommodate 
lepton flavour violation in the neutral lepton sector (i.e. neutrino oscillations), the contributions of 
these models to charged lepton flavour violating (cLFV) observables are typically extremely small. 
On the other hand, when such models - for example, the seesaw in its different realisations - are 
embedded within a larger framework, one can expect large contributions to cLFV observables, well 
within experimental reach. This is the case of super symmetric versions of the seesaw mechanism, 
which in addition to apport solutions to many theoretical and phenomenological issues, such as 
the hierarchy problem, gauge coupling unification and dark matter, can also account for neutrino 
data. 

However, these scenarios have several caveats, the most upsetting one being that they prove 
to be extremely hard to test, and thus can be neither confirmed nor excluded. This stems from 
the fact that in order to have sufficiently large Yukawa couplings (as required to account for large 
cLFV branching ratios), the typical scale of the extra particles (such as right-handed neutrinos, 
scalar or fermionic isospin triplets) is in general very high, potentially close to the gauge coupling 
unification scale. 

This can be avoided if one simultaneously succeeds in having TeV-scale mediators, while 
preserving the possibility of large Yukawa couplings. From an effective theory point of view this 
is equivalent to the decoupling of the coefficients associated to the dimension-five (at the origin 
of neutrino masses) and dimension-six operators (for instance, four-fermion operators): in other 
words, decoupling the smallness of the light neutrino masses from the flavour violation sources. 
For instance, this is possible in the case of the type-II seesaw (and its super symmetric (SUSY) 
realisations), as well as in the case of the so-called "inverse seesaw" (and the SUSY inverse seesaw). 

The inverse seesaw [8] constitutes a very appealing alternative to the " standard" seesaw reali- 
sations, and has recently been the subject of several dedicated studies [9H2T]. The inverse seesaw 
can be embedded in the Minimal Super symmetric extension of the SM (MSSM) by the addition 
of two extra gauge singlet superfields, with opposite lepton numbers (+1 and —1). As extensively 
discussed in [91-121] . in this framework one can in principle have large neutrino Yukawa couplings 
(Y v ~ 0(1)) compatible with a seesaw scale, M seesaw , close to the electroweak one (and thus 
within LHC reach), implying that there will be a significant enhancement of cLFV observables. 
This has fuelled a number of studies focusing on the potentialities of the inverse seesaw regard- 
ing cLFV and other phenomenological issues [9T (T2l[T4Hl7ll2TH2"4] . As recently discussed [25], the 
contributions of the comparatively light right-handed sneutrinos can enhance the Higgs-mediated 
penguin diagrams, leading to an augmentation of some observables - for instance Br(r — > 3/x) - by 
as much as two orders of magnitude, and have a non-negligible impact on Higgs-mediated leptonic 
B-meson decays and Higgs flavour violating decays. 
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These unique features - when compared to other SUSY seesaw realisations - open the door to 
rich phenomenological signatures, that can be potentially tested in the near future. The different 
high- and low-energy phenomenological implications of this class of models constitute the starting 
point to unveil the underlying mechanism of lepton mixing. 

There are currently a large number of facilities |26H34j . dedicated to the search of processes 
such as rare radiative decays, 3-body decays and muon-electron conversion in nuclei. Likewise, 
rare leptonic and semi-leptonic meson decays also offer a rich testing ground to experimentally 
probe cLFV. 

These (low-energy) searches are complementary to the LHC which, in addition to directly 
searching for new physics states, also allow to study numerous signals of cLFV at high-energy, 
typically in association with neutralino-slepton decay chains. In order to disentangle the under- 
lying model of lepton flavour violation, one relies on numerous strategies based on the interplay 
of low- and high-energy cLFV observables (see for example |35H38j ) . However, there are other 
avenues that can be explored in this quest to disentangle the underlying mechanism of neutrino 
mass generation, at the origin of lepton flavour violation: this approach is based upon exploring 
the correlation (or lack thereof) between different, unrelated, low-energy cLFV observables. The 
distinctive features of the underlying model will be manifest in the nature and specific hierarchy 
of the different contributions. For instance, in SUSY models where 7-penguins provide the domi- 
nant contribution to radiative and 3-body cLFV decays, one expects a strict correlation between 
Br(^ — > ej) and CR(/x — e, N). This is the case of constrained Minimal Supersymmetry Standard 
Model (cMSSM) based scenarios where additional lepton flavour violating sources have been intro- 
duced. Deviations from strict universality (as is the case of non-universal Higgs masses, NUHM), 
where for example Higgs-mediated penguins might play a significant role in fi — e conversion, break 
this strict correlation [39]. 

Independently of the specific mechanism of SUSY breaking, specific contributions to cLFV 
observables manifest a peculiar behaviour. A very interesting case is that of the Z-penguin: in 
the MSSM, the contributions of the Z-penguin to cLFV observables (such as ti — > 3£j and fj, — e 
conversion in nuclei) are suppressed by a subtle cancellation between the different terms in the 
amplitude. However, it has recently been noticed [30] that in models where new couplings are 
present or where the particle content is larger than that of the MSSM, the cancellation no longer 
holds, and the Z-penguin contributions can in fact provide the dominant contributions to cLFV 
processes such as 3-body decays, i{ — > 3£j and jjl — e conversion in heavy nuclei. Models with 
additional couplings (as is the case of trilinear R-parity violating super symmetric models) were 
discussed in 001131]. 

Contrary to the "standard SUSY seesaw", in the inverse SUSY seesaw the new states (in 
addition to those of the MSSM) do not decouple: indeed, the sterile states can be as light as 
to lie in the sub-GeV scale. In the framework of the inverse SUSY seesaw, one can see that 
the Z-penguins will also provide sizeable, if not dominant, contributions to a number of cLFV 
observables. 

In this work, we consider a realisation of the inverse seesaw, embedding it into an otherwise 
lepton flavour conserving supersymmetric extension of the SM, the cMSSM. We conduct a detailed 
study of the impact that enhanced Z-penguins might have on a large number of low-energy cLFV 
observables, in particular ^ — e conversion and /i — > 3e decay which in addition to being greatly 
enhanced by the Z-penguin, also have the best experimental prospects concerning the expected 
future sensitivities. Our results reveal that for vast regions of the parameter space, many cLFV 
observables are indeed boosted by the unsuppressed Z-penguin contribution and are within reach 
of present and future experiments. 
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Moreover, our analysis reveals that, similarly to what occurs in the MSSM, flavour chang- 
ing Higgs boson decays (where the Higgs boson contributions do not decouple with increasing 
supersymmetric masses |42H44| ) . the Z-penguin contributions to the LFV observables are not 
suppressed by a large SUSY scale. 

The paper is organised as follows: in Section [21 we will define the model, providing a brief 
overview on the implementation of the inverse seesaw in the MSSM. In Section [3J we discuss the 
enhancement of the Z-boson mediated contributions to low-energy cLFV observables. We also 
derive an analytical approximation for the Z — li — £j effective vertex. In Section HI we derive 
analytical expressions for several cLFV observables in the case where Z-boson penguin is the 
dominant contribution. In Section [5j we detail the corresponding numerical study, collect the 
relevant numerical results and discuss the results as well as the decoupling regime. Our final 
remarks are given in Section [6J 



2 Inverse seesaw mechanism in the MSSM 

The inverse seesaw model consists of a gauge singlet extension of the MSSM. Three pairs of singlet 
superfields, Df and X, (i = 1, 2, 3)0 with lepton numbers assigned to be —1 and +1, respectively, 
are added to the superfield content. The SUSY inverse seesaw model is defined by the following 
superpotential 



W = e ab 



jb , , fja fjb 



+ Y?9fL°H b u - vH a d H b u \ + M^Xj + -//.v. A\-.Y; , (2.1) 

where i,j = 1,2,3 are generation indices. In the above, and H u are the down- and up-type 
Higgs superfields, Li denotes the SU(2) doublet lepton superfields. The "Dirac"-type right-handed 
neutrino mass term M^. conserves lepton number, while the "Majorana" mass term Hx^ violates 
it by two units. In view of this the total lepton number L is no longer conserved; notice however 
that in this formulation (— 1) L remains a good quantum number. Since Mr^ conserves lepton 
number, in the limit ~* 0> lepton number conservation can be restored. In this study we 
consider a general framework with three generations of v c and X; we nevertheless recall that 
neutrino data can be successfully accommodated with only one generation of v c and X |17j . 
The soft SUSY breaking Lagrangian is given by 

-£ soft = -£^| SM + £fm§c m + xjm'X, + (^Y^^Zffl* + B^Mn^fXj 
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+ l Bi l i x ^x ij X i X j +h.c), (2.2) 

where £^| SM collects the soft SUSY breaking terms of the MSSM. B^ r and B% are the new 
parameters involving the scalar partners of the sterile neutrino states (notice that while the former 
conserves lepton number, the latter gives rise to a lepton number violating AL = 2 term). As- 
suming a flavour-blind mechanism for SUSY breaking, we consider universal boundary conditions 
for the soft SUSY breaking parameters at some very high energy scale (e.g. the gauge coupling 
unification scale ~ 10 16 GeV), 

= m , Mg aU gi no = M 1/2 ,A { = A I, B^ x = B Mr = B 1. (2.3) 



1 We use the notation: v c = v* R . 
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From Eq. (|2.ip one can verify that the two singlets vf and Xi are differently treated in the 
superpotential, so that, while a AL = 2 Majorana mass term is present for Xi (fiXi-XiXj), no 
[i v <}pfO^ term is included in W. Although the latter term can indeed be present in a superpotential, 
where (— 1) L is a good quantum number, we assume here for simplicity fi v ?. = 0. We notice that 
it is the magnitude of fix (and not that of fi u c) which controls the size of the light neutrino 
mass |13pi8| , and that the absence of the mass term fi v c does enhance the symmetry of the model 
(a non- vanishing, but small value of fi v c. does not have any effect on the qualitative features of 
the model). Here, fix will be considered as an effective parameter, no assumption being made on 
its origin. Such a parameter could be understood either dynamically or in a SUSY Grand Unified 
Theory framework [HllIIllIHi 20j . Furthermore fx v c <c fix can also be realised in some extended 
frameworks |13j . 



2.1 Neutrino masses 

We consider a general framework with three generations of sterile neutrinos uf and Xi. Conse- 
quently, one has the following symmetric (9 x 9) mass matrix Ai in the basis {v,v c ,X}, 




M = m D M R , (2.4) 



Where, mo = ^*Y w v u and Mr, fix are (3 x 3) matrices in family space. Assuming mp, fix "C Mr, 
the diagonalization leads to an effective Majorana mass matrix for the light neutrinos 



m v = m T D Ml X fi x M^m D = "-^Y'i (M^)" 1 fi x M^Y v . (2.5) 

As mentioned earlier, one of the advantages of the inverse seesaw mechanism is that the smallness 
of the light neutrino masses is directly controlled by the size of fix- Hence the lepton number 
conserving mass parameters ttid and Mr can easily accommodate large (natural) Yukawa cou- 
plings (Y u ~ 0(1)) and a right-handed neutrino mass scale around the TeV, see Eq. (|2.5p . In 
turn, this allows to have sizable contributions to cLFV observables, contrary to what occurs in 
the framework of, for example, a type-I seesaw. 

In analogy to a type-I seesaw, one can define an effective right-handed neutrino mass term M 
such that 

M" 1 = (Mr)' 1 fix Mr 1 . (2.6) 

With this definition, the light neutrino mass matrix can be cast in a way strongly resembling a 
standard (type-I) seesaw equation 

m v = ^yjM- 1 ^. (2.7) 

This effective light neutrino mass matrix (m u ) can be diagonalized as 

^PMNS m ^PMNS = diag mi . (2.8) 

Using the above equations, one can express the neutrino Yukawa couplings (Y u = ^mo) as in [9] 
(in analogy to to the Casas- Ibarra parameterization |46| in standard seesaw), 

y v = ¥-Vm Ry/m c4 MNS , (2.9) 



v 
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where V M is the diagonal matrbQ with the square roots of the eigenvalues of M and R is a 3 x 3 
orthogonal matrix, parameterized by 3 complex angles, which encodes additional mixings. In our 
study of the different cLFV observables, we assume specific regimes for R. 

Without loss of generality, we choose a basis where Mr is diagonal at the SUSY scale, i.e., 

M Rij = ding M Rii . (2.10) 

In addition, in the numerical evaluation, we shall also assume nx^ to be diagonal, a simplifying 
assumption motivated by the fact that cLFV observables depend only indirectly on [ix^ , as already 
discussed in the introduction. We will further explain this issue in Section 15.11 In the subsequent 
analysis, we assume Mr, fix to be free parameters, determined in order to comply with neutrino 
data. 

Concerning neutrino oscillation data, there has been a very intense experimental activity re- 
lated to the #i3 mixing angle, with very recent results from Double-Chooz PQ, T2K [2j, MINOS [3], 
Daya-Bay [I] and RENO [S]. We use in our analysis the (best-fit) results of [37], supplemented 
with the Daya-Bay measurement for #13 [4]. Therefore, we consider the following set of parame- 
ters, namely the two mass-squared differences and the three mixing angles, given below with 3a 
range [4ll47|. 

sin 2 0i2 = 0.27 - 0.36 , 

sin 2 <9 23 = 0.39 - 0.64 , 

sin 2 2#i3 = 0.092 ± 0.06 , 

Am|i = 7.09 - 8.19 x 10" 5 eV 2 , 

lAm^l = 2.18 - 2.73 x 10~ 3 eV 2 . (2.11) 

We also assume normal hierarchy (Am^ > 0). Updated global fits to all available experimental 
data have also appeared recently [SUTJ. However, these do not have a significant effect on our 
numerical results. 

3 Z-boson mediated cLFV 

In this section we shall discuss the enhancement of the Z-boson mediated contributions, presenting 
approximate results for the Z—£i — £j effective vertex. We then proceed to discuss how, in different 
extensions of the MSSM, the Z-boson mediated contributions can significantly enhance several 
cLFV observables. The latter effect, which is absent in the MSSM, can have a strong impact on 
the theoretical predictions for cLFV rates, potentially leading to a very different phenomenology 
since Br(£j —> ij-y) will no longer be the most constraining observables. 

3.1 Enhancing cLFV with the Z-penguin 

In the MSSM and in many of its extensions, photon penguins provide the dominant contributions 
to 3-body cLFV decays ti — > 31 j [481149], The only exception arises in the large tan/3 (and low 
tua) regime, where large Higgs contributions are expected |5offi However, it has recently been 
shown that many simple extensions of the (s)leptonic sector can lead to large enhancements of the 

2 In the subsequent sections we always consider scenarios where M = diag(M, M, M). Therefore, for the sake of 
brevity, we use the simple notation M both for the matrix and its eigenvalues. 

3 A similar observation in the context of the inverse seesaw was recently made in [25] - 
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.Z-boson contributions 00], potentially leading to scenarios where the latter provide the dominant 
contribution to Br(£j — > 3£j) (> Br(^j — > ^7)). This can be understood from simple mass scaling 
considerations: let us consider the chargino-sneutrino 1-loop diagrams leading to ii — > 3£j. The 
photon-penguin contribution can be written as 



A^ R = t-^O l a > R s(x*) , (3.1) 



16vr 2 



2 

m~. 



whereas the Z-contributions read 



Fx = ~ 2 2 . 2fl zO L F ft(x 2 ) , (3.2) 

lo-K z g z sin Owf^z 

with X = {LL, LR, RL, RR}. In the above Oy' R denote combinations of rotation matrices and 
coupling constants and s(x 2 ) and t(x 2 ) represent the Passarino-Veltman loop functions which 
depend on x 2 = m 2 _/m 2 (see [19]). Notice that the only mass scale involved in the A form 
factors is mgusY (the photon being massless). On the other hand, the mass scale in the Fx form 
factors is set, in this case, by the Z-boson mass (mz)- Therefore, we conclude that A ~ ""T-gusy 
and F ~ m ~z 2 - Since m 2 z <C ^lusY' the Z-penguin can, in principle, dominate over the photon 
penguin. Assuming that all loop functions, mixing matrices and coupling constants are of the 
same order, one can estimate 

F m 2 

T ~ 2 - 2 V a Y 2 ~ 500 for ^usy = 300 GeV. (3.3) 
-r± cj sin u'Y/Vfiy 

Moreover, the ii — > 3ij decay width depends on F 2 and A 2 , and thus the above ratio becomes even 
more pronounced!^. However, a subtle cancellation between the different diagrams contributing 
to the leading ^-contribution |40| implies that in the MSSM, the photon penguin is found to be 
(numerically) dominant |48|[49j. To understand this, notice that the dominant contribution to 
ii — > 3£j comes from diagrams where the leptons in the external legs are left-handed (the other 
cases are suppressed by the Yukawa couplings of the charged leptons). This is given by the form 
factor Fll, usually written as 



where zjp = — ■£-{— \ + s w) 1S ^ ne Z — ii — ij tree-level coupling (i = j) and Fl is the Z — ii — ij 1- 
loop effective vertex, with i ^ j, and with cw = cos 0y/ and sy/ = sin 6\y. Fl receives contributions 
from different 1-loop diagrams, and here we focus on the chargino-sneutrino loop contribution. 
Expanding in the chargino mixing angle, , one can write (see Figfl]) 

F L = Ff + \0\ ± F^ + ... . (3.5) 

Notice that there is no term in the expansion involving at the leading order, nor at the 1st 
order, since there is no — vl ~ coupling. For this reason, only the wino contributes at the 
zeroth order in 0^±. F^ can be written aa^| 

*f% = F ® = "IS? ( M ^ ve + M ^ + M » 2 ) ' (3 ' 6) 



4 From these considerations one can also conclude that in a very light SUSY scenario, with msuSY ~ mz, photon 
and ^-contributions to ii — > 3ij are of the same order. This observation has been confirmed numerically, see section 

5 For the sake of clarity, we omit flavour indices in the expression for Fl (and in the expansion coefficients). 
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(a) Oth order term: (b) 2nd order term: 

Figure 1: Diagrammatic representation of the Fl expansion in the chargino mixing angle, 0^±. 

with 

M^ ve = ~g 2 (gc w -g's w )Z v i *Z^f^ ve , (3.7) 

M% = g 3 c w ZrZ^f^, (3.8) 
1 
2' 

where a sum over the index x is implicit. The terms in the sum come from different types 
of diagrams: wave function diagrams (M wavc ) and penguins with the Z-boson attached to the 
chargino line (M p i) or to the sneutrino line (M p 2). Moreover, Zy is a 3 x 3 unitary ma- 
trix that diagonalizes the mass matrix of the sneutrinos. Here /^ ave = —Bi(m 2 ± ,m 2 ^),f^ = 

\C {ml x ,m 2 ^,m 2 _±) - ml±Co(ml x , mi±, m|±), f% 2 = ^(m^m^rr^); for the exact defini- 
tions of these loop functions, see [39]. The sum in Eq. (|3.6p exactly vanishes, as can be verified 



M% = --g 2 (gc w + g's w )Z v i *Z v j f^, (3.9) 



by grouping the different terms 

if ) = ^(a cwZrz^Xf + g's w zrz^xi) (3.10) 

with Xf = /^ ave - 2/£ + /p 2 ,X| = f£ 2 - /* ave . Using the exact expressions for the loop 
functions one finds that the masses cancel out and that these combinations become inde- 
pendent of x = m-f-fmp' Xf = X\ — — | and X 2 = X 2 = —\, Vx. Therefore, one is left with 

if'* oc Zy *Zy = (ZyZy\ , which vanishes for i 7^ j due to unitarity of the Zy matrb<[^|. 
In conclusion, the first non-vanishing term in the expansion appears at 2nd order in the chargino 
mixing angle, which naturally suppresses the Z-mediated contributions. This is the reason why 
the photon contributions turn out to be dominant in the MSSM. 

However, there are many cases where the cancellation of the zeroth order term in the expansion 
no longer holds, as discussed in [40j . where numerical examples were given, including the inverse 
SUSY seesaw. The introduction of new interactions for the (s)leptons, in particular, the one 
involving the Yukawa couplings Y u , modifies the previous conclusion for the MSSM: owing to 
the — ur — £l coupling, higgsino contributions for if ^ preclude the previously discussed 

cancellation. In fact, there is a non-zero if ^ contribution that enhances Fl by a huge factor 
leading to large Z-penguin mediated contributions to cLFV processes. 

It is worth mentioning that the discussed cancellation only happens for wino diagrams and not 
for the higgsino ones. Therefore, Fr has already a finite zeroth order contribution in the MSSM. 
However, this is very small because of the tiny charged lepton Yukawa couplings. 



3 We notice that a similar behaviour was found in [5T] in association to the B — > X s £ + £ decay. 
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Finally, although the previous discussion has been focused on 1^ — > 3£j, the same enhancement 
in the Z — t{ — lj effective vertex also affects other observables which are mediated by Z-boson 
exchange. This is the case for — e conversion in nuclei [39J and r — > P°£i, where P° is a 
pseudoscalar meson [52] . 



3.2 Approximate expression for the 1-loop Z — — £j effective vertex 

This section is devoted to deriving an analytical approximation for the 1-loop Z — ii — £j effective 

vertex, Fl — F^ , in the framework of the SUSY inverse seesaw. 

Following the previous discussion, we will define this vertex with the lepton li as incoming 
particle. The expression for F^ obtained from the chargino-sneutrino loopd3 can be decomposed 
as 

r L — 167r 2 \^ L,wino ' L,higgsino J J \°- LL ) 

where F^- mo and -F^ ^ ino are the wino and higgsino contributions, respectively. The latter 
terms can be rewritten as 

-2 



*£Lo = ~Y Z v* Z v (acwY? + g'swYf) (3.12) 

. r / _ i \ 1 , , „ , „ ~ _ _ .i 

(3.13) 



L.higgsino 4 v v 



1 1 1 jf q 7 y,*'+3* 72,2+3^ / 2 ^2 3 



where Zy is now a 9 x 9 unitary matrix and sums over x,y = 1,...,9 and z,z' = 1,2,3 are 
implicit. We have also defined the following combinations of loop functions: 

Ff = -C (m^,m| ± ,m|±) + 2(7o(m^,m| ± ,m| ± )-fi 1 (mL,m?J (3.14) 
Yjt = ^S xy C„{m^,ml,m%) + Biim^miJ (3.15) 

and 

3 

S xy = Y / Zv* z v- ( 3 - 16 ) 

fe=i 

Equations (|3. 12[) and f|3. 13|) are exact and do not involve any approximation nor assumption 
on the sneutrino mixing pattern. However, in order to render these expressions more transparent, 
let us now consider the following limit, 

c _ / $xy forx,y < 3 . s 

forx,y>3 

This assumption actually corresponds to the MSSM, since in this case Zy is a 3 x 3 unitary 
matrix. It also provides a very good approximation in the inverse seesaw when the Yukawa 
couplings are relatively small. In this framework^ the sneutrino mixing matrix Zy , written in the 



7 The analogous neutralino-slepton loop diagrams provide extremely small contributions since they involve 
charged lepton Yukawa couplings and small mass insertions on the slepton propagators. 

8 As will be shown in our numerical analysis, regimes of small Yukawa couplings are actually favoured by current 
bounds on cLFV observables. For instance, /i — e conversion in gold leads to {YjY„) 12 < 3 x 10" 4 . 



9 



basis v x = (^1,2,35^123)^1,2,3)) is approximately given by 

(3.18) 



h 











where I3 is the 3x3 identity matrix and Zy is a 6 x 6 unitary matrix that diagonalizes the singlet 
sector. The negligible left-right mixing in the sneutrino sector thus justifies the approximation of 

Eq. una). 

In both cases, MSSM and the inverse seesaw with u mixings as defined in Eq. (|3.18p . one can 
further simplify Eq. f)3. 12j) by means of Eq. (|3.17j) . It is straightforward to verify that Yf reduces 
to Y* = X\ = — |, whereas Y% reduces to Kf = X2 = — \- Therefore, using unitarity of the Zy 
matrix, Zy* Zy = 5ij, the wino contribution simplifies to 

F LMno = Y Si ^cw + 9'sw) = ^<%(1 + 24) , (3.19) 

which vanishes in the case of flavour violating transitions (i ^ j). This is in perfect agreement 
with the Z-penguin cancellation in the MSSM, discussed in the previous subsection. 

Let us now consider the higgsino contributions F^l^ no * P resen t in the MSSM. In this 
case, the results also simplify under the assumption of Eq. (|3.18p . which allows to cancel the last 
term in Eq. ()3.13p . Therefore, ^°higg S i no can De written as 

*2L- = 4 (tf*), ,(*- 5) • <3 - 20) 

The above equation corresponds to the approximate LFV Z — li — £j 1-loop effective vertex. 

We have explicitly checked that this formula does indeed reproduce the full numerical results to 
a very good approximation. Further refinements can be obtained by including other relevant (but 
sub-dominant) 1-loop diagrams. In particular, we found non-negligible sub-leading contributions 
arising from 1-loop diagrams involving charged Higgs and neutrinos. We also note that 
and the higher order terms in the expansion receive new contributions. However, the latter are 
(numerically) negligible when compared to 

Finally, it is worth noticing that the effective vertex exhibits a crucial property: at 
leading order there is no dependence on any supersymmetric mass (however higher order terms in 
the expansion do indeed decrease for increasing SUSY masses). 

There are well-known results in general quantum field theory (for both SUSY and non SUSY 
cases) regarding the decoupling theorem [53^154] . which seems not to apply in this case. Z- 
mediated processes exhibit a non-decoupling behaviour and large supersymmetric masses do not 
suppress the charged lepton flavour violating signatures induced by Z-boson exchange. This non- 
decoupling behaviour has also been found in flavour changing Higgs boson decays in the MSSM, 
both to hadronic |42y43| and leptonic final states 



4 Lepton flavour violating observables 

4.1 Current experimental situation and future prospects 

The search for cLFV is a very active field with either dedicated experiments like MEG [55] or others 
with a broader program like B factories |56| . In this paper, we focus on leptonic observables, which 
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can be classified as radiative decays, e.g. \x — > &y, 3-body decays, e.g. r — > 3/i, and neutrinoless 
conversion in muonic atoms, e.g. fx, Ti — > e, Ti. 

The experiments looking for radiative decays are quite different depending on the lepton in the 
initial state. If it is a muon, the only decay is /i — > which is studied by dedicated experiments 
such as MEG |55] , This collaboration has plans for an upgrade that would improve the sensitivity 
to O(10~ 14 ). Radiative r decays are studied at B factories, which are also r factories, since 
the production cross-sections are very close at the T(4s) resonance. The current upper limits 
on Br(r — > /X7) and Br(r — > ej) are given by the BaBar experiment, together with expected 
sensitivities at the future generation of B factories, e.g. Belle II and SuperB [56] , 

For the same reasons 3-body decays of the r lepton are also usually searched for at B facto- 
ries. The current upper limits come from the Belle experiment [571158] because of its larger data 
sample compared to BaBar. Since these observables are currently not limited by the background, 
significant improvements are expected at Belle II and SuperB [56]. The decay \i — > 3e has been 
investigated by the SINDRUM experiment |28j and, if approved, a future experiment named Mu3e 
at PSI could reach a sensitivity of 10~ 15 (after upgrades 10 -16 ) [59] , 

Neutrinoless conversion in muonic atoms has also been studied for different nuclei by the 
SINDRUM II collaboration [601161] which has set the current upper limits. In the future, the 
sensitivity is expected to be greatly improved by different project^]. For convenience, in our 
numerical study we will consider a future sensitivity in the range HT 16 - 10~ 18 . 

In addition to these low-energy observables, interesting phenomena are expected to be observed 
at colliders. One can have sizable widths for processes like X2 ~* Xi^i £f 1 flavoured slepton mass 
splittings (especially between the first and second generation of left-handed sleptons) and finally 
the appearance of new edges in same-flavour dilepton mass distributions. Assuming a unique 
source of LFV (neutrino mass generation) , the interplay of low- and high-energy LFV observables 
can strengthen or disfavour the underlying model of new physics. Illustrative examples of the 
potential of this interplay can be found for instance in [35H38] . 



4.2 \i — e conversion in nuclei 

This process is particularly sensitive to the enhancement of the Z-penguin contribution. The 
conversion rate, relative to the the muon capture rate, can be expressed as [39] 



CR(/i — e, Nucleus) 




Z and iV are the number of protons and neutrons in the nucleus and Z e g is the effective atomic 
charge [67]. Similarly, F p is the nuclear matrix element and r cap t represents the total muon 
capture rate. Gf is the Fermi constant, a is the fine structure constant, p e and E e ( ~ in the 
numerical evaluation) are the momentum and energy of the electron and m M is the muon mass. 

9 Mu2e [621163] is a future experiment at Fermilab with expected sensitivities of respectively 10 -17 (phase I) and 
10 -18 (phase II with Project X). On the other hand, the first experiment that could be built at J-PARC is DeeMe [64] 
with an expected sensitivity of 2 x 1CT 14 in 2015. Then COMET [65] and PRISM/PRIME [66] would come with 
sensitivities of 10" 15 (COMET Phase I, 2017), 10~ 17 (COMET phase II, 2021) and 10" 18 (PRISM/PRIME) for a 
titanium nucleus. 
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In the above, g XK anc ^ 9xk ( W1 th X = L,R and K = S, V) represent the relevant isoscalar, and 
isovector couplings respectively, and are given by 



(<?,«■) 

K 



q=u,d,s 

2 ^ 

q=u,d,s 



(0) _ 1 / r (q,p) r 

9XK ~ 2 2-> \9XK(q)^K + 9XK{q)^ 

riii'P) „ Mi' 71 }} 

9XK{q)^K ~9XK{q) (j K J 



(4.2) 



The numerical values of Gk can be found in |39j. 

Similar to other observables involving four fermions, the [i — e conversion rate receives con- 
tributions from 7-, Z- and Higgs-penguins as well as box diagrams. The corresponding couplings 
are, 



7 Z B 

9LV(q) = 9 L V{q) + 9LV(q) + 9LV{q) > 

H B 
9LS{q) = 9lS(q) + 9lV(q) 



(4.3) 



dominant contribution 10 !, we focus on gfy^ q y We thus have 



where gl x(q) , Q Z LX ( q y sfs( g )> g LX( g ) ( with X = V,S) represent the couplings of the photon, Z, 
H and box diagrams, respectively. Again, considering that the Z-boson is at the origin of the 



(4.4) 
(4.5) 



9LV(q) = 9 L V(q) 



G F 2 



rrir 



9\Q) — ' 



and gRv( q ) = 9hV(q) \l^~>r- The Z — q — q couplings Z q R ) can be written as 

z\ 



J L 
J R 



CW 

Wq b W> 

CW 



q a W\ i 



(4.6) 
(4.7) 



with sw = sin#w and cw = cos 6w- Then, considering that Fr <C -Fl, we have 



Ve E e mlG l F o?Z^Fl 



8tt 2 Z 

,3 n2 ™3 7 4 



{Z + N)gfl + {Z-N)g 



Pe E e mtGj,at i Z* s F 1 



\ir 2 Z 



(Z + N)g L % + (Z- N)g™ 



CR([i — e, Nucleus) 



where 



q=u,d,s 

Finally, substituting Fl and denoting the hadronic coefficient 
CHad, we have 

Pe E e mlG 2 F a 3 Z^F 2 g 2 ^ (c w 



(i) 
LV 



capt 



(1) 



9^ = 1 E s(*)(<# rt ±Gfr ) 



2 Ft 



(4- 



(Z + N)g® +{Z- N)gM 



CR(/i — e, Nucleus) 



2 17 ^Zc w 



X CHad" 



' )2 W 



1 



12 rv 



(4- 



<^ ' j- capt 

This approximate formula is valid for scenarios where Z-boson provides the dominant contribution 
to the conversion rate. 



We stress that in our numerical analysis, we will take into account all contributions to the CR. We refer the 
reader to [35] for the corresponding formulae. 



12 



4.3 Br(£i -> 3^) 

As we discussed in Section the Z-penguin can provide the dominant contribution to the decay 
width of t{ — > 3£j in the inverse seesaw extension of MSSM. In the limit of Z-penguin dominance, 
the decay width can be cast as (SHIP 



Uj) ~ 5i^ m 'J (i? ^ + Frr + f2lr + Frl) 



(4.10) 



_ e _„5 z 772 



The complete expression for Fll can be derived using Eqs. (|3.4p . ([3.1ip . (|3.19p and (|3.20p . 



Fll 



(2 s w) 



1287r 2 s^c^m| 

„2 \2 



-2 S VK^ 



l2S-K 2 s 2 w c 2 w m? z 



g 2 5 ij (l + 2c 2 w ) + (Y^Y l/ )lc 2 w 



so that the decay width for r(£, — >• 3£j) is given by 

— >■ 3£j) = 22 2 7 4 4 ^ m fi [YyYv 



(4.11) 
(4.12) 



(4.13) 



In this limit one recovers a correlation between CR(/i — e, Nucleus) and Br(/j — > 3e), which we 
shall briefly mention in our numerical results. 



4.4 Br(r P ^) 

r decays into a light charged lepton and a pseudoscalar meso P° = 7T°, 7/, r/, can also be signif- 
icantly enhanced in this framework, owing to Z-boson mediation. The corresponding branching 
ratio can be written as [52] 



Br(r -> P° 



4vr m?r T 2 ^ 1 1 ' 



(4.14) 



where A(x, y, z) = {x + y — z) 2 — Axy and T T is the total decay width of the r lepton. The averaged 
squared amplitude summed over initial and final states is given by 

I \T\ 2 = [2m^.m T (apaf* - b k P bfA + (m 2 T + m| - m|) fea^* + 

fc,m 



with k,m = Z,A°. Focusing on the Z-boson contributions, we have 

9 Pr C(P) 



(4.15) 



dp 



bp 



2 cos 2 

g P. g(P) 
2 cos #w 2 m \ 



(m T - me.) (Pl + Pr) 



(4.16) 



1 We do not consider final states with a if as this would involve additional suppression due to the double penguin. 
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Figure 2: Feynman diagram contributing to Br(x2 — > Xi^j)- The white white circle represents a 
LFV mass insertion in the charged slepton propagator. 

where, Fr = Fl\l++r, which is in general small. In the above, F w ~ 92.4 MeV is the pion decay 
constant while the functions C{P) have been defined in [52] . In the limit Fr <C Fl, we have 

I E \ T \ 2 = n J g2j f C( f )2 [-8m?m|. + (m 2 . + ml - m P )(m 2 T + mj.)) Ff . (4.17) 
Finally, using the expression for Fl, one finds 




(4.18) 



Although this approximate formula corresponds to leading order estimates, in the numerical eval- 
uation, we take into account all contributions from the Z- and Higgs bosons. 

4.5 Br(l« -> Iff) 

For completeness, we discuss the radiative decays — > if. These observables, which have been 
studied in great detail in [9"j l48l lS8]. are not sensitive to the Z-mediation. It is however known 
that, for small values of Mr, the branching fraction can easily reach the experimental sensitivity, 
even in the absence of supersymmetric contributions In the following section, considering both 
heavy singlet neutrino as well as supersymmetric contributions, we will explore the part of the 
parameter space where Br(^j — > if) is within experimental sensitivity reach. 

4.6 Collider observables: Br(%2 — > Xi^dj) an d Am f - 

So far, we have only addressed low-energy observables. However, cLFV collider observables such 
as Br(%2 - > Xi^i^j) anci the slepton mass splittings Am f - (m £ ~ — mg) are also relevant in scenarios 
that provide a strong connection between low- and high-energy lepton flavour violation |35H38j . 

However, in the model under investigation, one expects that these high-energy observables 
will be of little relevance when compared to the low-energy cLFV observables. For example, the 
decay \2 ~* Xi^i^j 1S induced by diagrams like the one shown in Figure [2j As we will verify in the 
following section, the experimental limits on the low-energy observables require small neutrino 
Yukawa couplings (typically of order ~ 10 -3 — 10 -2 ), which in turn leads to low ^2 ~~ ^ Xi^dj rates 
and small mass splittings. 

For completeness, we have also studied LFV decays such, as Z — > iiij and h — > iiij. However, 
as explained below, no enhancement was found. 
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5 Numerical results 



In this section we present our numerical results. We begin by introducing the basic setup for our 
computation and show some results concerning the relative importance of the different contribu- 
tions to the cLFV observables. We then discuss our results for the rates of cLFV observables, 
emphasising the most relevant features of this model. Finally, we provide a brief summary of the 
main results, which we illustrate via some representative benchmark points. 

5.1 Setup 

Our numerical results have been obtained with a SPheno [691170] code generated with the Math- 
ematica package SARAH [71H73] . The computation of the LFV observables is based on the results 
presented in [39,49,52 , which we extended to the inverse seesaw case. 

In what concerns the supersymmetric parameters, the model is characterised by 



Due to their little influence on the relevant observables, we fix sign(/x) (+) and Bq(= 0) in our 
scanJ^l. We have explicitly checked that only /i — > ej depends on the Bo parameter. This can be 
understood from its effect on the sneutrino spectrum, which in turn strongly affects Br^ — > ej). 
We find that in some cases one can enhance Br(^ —> ej) by as much as one order of magnitude if 
one uses Bq to fine-tune the superparticle masses. The other observables are nearly independent 



Concerning the remaining cMSSM-like parameters, mo, M1/2, Aq and tan/3, we will consider 
some specific scenarios. We nevertheless recall that the Z-penguin, dominant in most of the 
parameter space, has a very mild dependence on the cMSSM parameters, so that our conclusions 
for the Z-boson mediated processes are quite general. On the other hand, for the radiative 
li — > £jj, the dependence on the SUSY parameters cannot be neglected. These observables have 
been studied in SUSY and non-SUSY scenarios [9l[TT l l?4" l lT5] and we include them in our analysis 
for completeness. As we will see, only for very low values of Mr or when the superparticles in the 
loop are light, they become the most constraining observables. 

As explained in Section 12.11 we will work in the basis where Mr is a diagonal matrix and 
compute the resulting Yukawa couplings by means of the Casas-Ibarra parameterisation [46] . In 
this basis, both Y u and fix can be in general non-diagonal. In fact, neutrino mixing requires 
off-diagonal entries in at least one of these two matrices. In principle, one can consider two limits: 
(1) Y v limit: Y v contains off-diagonal entries and fix is diagonal, and (2) fix limit: Y u is diagonal 
and fix contains off-diagonal entries (any phenomenologically allowed scenario will be either Y v 
limit, fix limit or an intermediate case). However, the off-diagonal elements in fix have no impact 
on the phenomenology in the charged lepton sector, since this parameter does not appear in the 
coefficients of the operators relevant for the cLFV transitions. Therefore, in the fix limit (where 
Y v is diagonal), all cLFV observables vanish. Thus, we choose to work in the Y v limit, assuming 
fix = fix I; where I is the 3x3 identity matrix, which in turn maximises the lepton flavour 
violating effects allowed by this model. However, we stress that the model cannot be ruled out by 
the non-observation of cLFV processes, since one can always approach the fix limit to suppress 
the corresponding cLFV observables. 

Concerning the Mr matrix, we will consider two scenarios: degenerate and non-degenerate 
spectrum, for simplicity presenting our results for the first one. In any case, and as shown 

12 It is worth recalling that Bo only affects the sneutrino sector, fi and £? M are computed using the tadpole 
equations at the SUSY scale. 



m , M 1/2 , A) , tan /3, sign(^) , B . 




oiB Q . 
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below, the Z-boson mediated processes have very little dependence on the right-handed (s)neutrino 
spectrum. We further choose R = 1 (see Eq. (|2.9p ) in order to keep the discussion as simple as 
possible^!. 

Let us now discuss the relevant parameters for the cLFV observables. As explained above, 
fix has very little impact on the phenomenology, since it only affects LFV due to its relation to 
the size of Y v . On the other hand, Mr has a stronger impact, since it affects the masses of the 
singlet states, but, as already explained, the Z-penguins are only slightly sensitive to the spectrum. 
Therefore, the impact of Mr on cLFV is rather small, only via the size of the Yukawa couplings. 

The relevant quantities for the analysis are M ~ — a (as defined in Eq. (|2.6p ). controlling the 

i — i 

size of the Yukawa couplings Y u , and Mr, which plays a sub-dominant roleo Therefore, we will 
study the variation of the cLFV rates with respect to M, showing how different possibilities for 
Mr affect the numerical results. 

In this work, we do not take into account the LHC Higgs mass constraint on the allowed 
parameter space. Our conclusions for lepton flavour violating observables are expected to hold 
when the Higgs mass constraint is forced upon the cMSSM parameter space. Indeed, one can 
always find regions in the cMSSM parameter space that fulfill the recent results from the LHC 
[761177] and Tevatron [78], pointing towards a Higgs mass in the 125 GeV range: this can be 
achieved with large Aq (to increase the mixing in the stop sector) and a moderately large tan f3, 
see for example [79H81] . This conclusion can easily be generalized to the supersymmetric inverse 
seesaw. 

In the limit of small Yukawa couplings, the singlet sector contributions to the Higgs mass 
are negligible [82]. In fact, and as we will see below, limits on cLFV processes set important 
constraints on the Yukawa couplings, forcing them to be below Y v ~ 10~ 3 — 10~ 2 (unless one is 
working in the \xx limit defined above). Thus one simply recovers the cMSSM result. We will 
discuss a particular benchmark point where this is explicitly shown. 

5.2 Relative importance of the different contributions 

We start this section by considering the different contributions of the SM/SUSY particles to the 
cLFV observables. We display our results for 1^ — > 3£j and \x — e conversion in gold, although the 
conclusions drawn here can be easily extended to other cLFV observables. 

Figure [3] shows how the different relative contributions to Br(/z — > 3e) and \i — e conversion in 
79 7 Au depend on the SUSY scale, setting for simplicity mo = Mi/ 2 . The relative contributions, 
denoted defined as 

c t = <glL_ , (5.2) 

where Rj is the rate (branching ratio in case of fi — > 3e and conversion rate in case of \i — e 
conversion in Au) that would be obtained should the z-contribution to the process be the only 
one. These numerical results have been obtained with fixed values for Aq = —300 GeV, Bq = 0, 
tan/3 = 10 and sign(//) = +. Nevertheless, we have verified that a different set of parameters 
would lead to similar results. 

For Br(^i — > 3e) the different contributions include (7, Z, h, interference 7 — Z, interference 7 — 
h, interference Z — h). Note that h includes both Higgs and box contributions, grouped together 



13 In general, the limit R — 1 translates into a "conservative" limit for flavour violation: apart from possible 
cancellations, and for a fixed seesaw scale, this limit typically provides a lower bound on the generated cLFV rates. 

14 The exception to this rule is fi — ¥ e^f. as one lowers Mr, one enhances Br(/i — > ej); this is due to non-SUSY 
contributions, see Section [4.51 Such an effect is less pronounced in other observables. 
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500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 

mo = Mi/ 2 [GeV] mo = M 1 / 2 [GeV] 

(a) fi — > 3e (b) jj, — e conversion in 7g 7 Au 

Fi gure 3: Relative contributions to Br(/i — y 3e) (left-hand side) and fi — e conversion in ^g^Au 
(right-hand side) as a function of mj = Mjy 2 for a degenerate singlet spectrum with Mr = 10 
TeV and M = 10 11 GeV. The rest of the cMSSM parameters are set to A Q = -300 GeV, B Q = 0, 
tan/3 = 10 and sign(^) = +. On the left-hand side, solid lines represent individual contributions, 
7 (black), Z (blue) and h (red) whereas the dashed lines represent interference terms, 7 — Z 
(green), 7-/1 (purple) and Z — h (orange). Note that in this case h includes both Higgs and 
box contributions. On the right-hand side interference terms are not shown to make the results 
clearer. The individual contributions are 7 (black), Z (blue) and h (red) and boxes (purple). 



in [J9]. Figure 3(a) clearly shows that the Z-boson contribution to the process is the dominant one. 
Only for very low mo = M\/2 can one find competitive (even dominant) photon contributions, as 
expected from the theoretical arguments, see Section [3j However, the low mo region has already 
been excluded by direct collider searches [831153]. 

A very similar behaviour is found for fi — e conversion in Au, where the Z-boson contribution 
turns out to be the dominant one as well. Figure [3 (b ) | depicts the different individual contributions 
to the conversion rate, with i = (7, Z, h, boxes). In this case, interference terms have not been 
shown to simplify the plot and to render the main results more visible. Again, as occurred for 
Br(/i — > 3e), the photon contribution becomes important for light SUSY scenarios (lower values 
of mo = Mi/2)- The little dip in the boxes curve is due to small numerical instabilities in the 
loop function computation. The reduction of the photon contribution in fj, — e conversion, when 
compared to the Br(/i, — > 3e), is due to the smaller electric charges of the quarks compared to 
those of the leptons. 

Finally, we show in Figure |4] the absolute individual contributions to Br(/i — > 3e) as a function 
of mo = Mi/2- This plot, complementary to the one on the left-hand side of Figure [31 serves to 
illustrate the non-decoupling behaviour of the Z-boson contributions. As one increases mo = Mi / 2 , 
the resulting SUSY spectrum becomes heavier and the photon, Higgs and box contributions to 
Br(/i —7- 3e) clearly decrease. In contrast, the Z-penguin contribution remains constant in the 
heavy SUSY limit, in which case the 1-loop effective Z — [i — e coupling can be well described by 
the leading order term , see Eq. (|3.20p . 

In conclusion, in the framework of the inverse seesaw, should the Z-penguin be present, it will 
provide the dominant contribution to cLFV observables (except in the case of a very light SUSY 
spectrum). This result clearly justifies the approximations of Section UJ 
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500 1000 1500 2000 2500 3000 
m o = Mi/2 [GeVJ 



Figure 4: Absolute contributions to Br(/i — > 3e) as a function of mo = M\/2- The parameters and 
the color code are the same as in figure [3j 

5.3 /i — e conversion in nuclei 

After the discussion on the relative size of the different individual contributions, we proceed to 
study the different cLFV observables and begin with [i — e conversion rates in nuclei 15 !. 

We first address the degenerate singlet scenario, Mr = diag(M#, Mr, Mr). The more general 
case of a hierarchical spectrum will be briefly discussed in Section f5. 7\ although we should mention 
that the difference for this observable is small and does not affect our conclusions. 

Figure [5] shows the fi — e conversion rates for different nuclei 7g 7 Au (left) and IfTi (right), as a 
function of M for three different Mr values (100 GeV, 1 TeV, 10 TeV). Again, we take A = -300 
GeV, tan j3 = 10, sign(//) = + and Bq = 0, but mo and Mi / 2 are randomly taken in the range [0, 3 
TeV]. Note that there is a sharp correlation with M, hardly distorted by the changes in mo and 
Mi/2- For the red and black dots, this is a consequence of Z-boson dominance, as discussed in 
Section O The blue dots correspond to the limiting case of non-SUSY photon-penguin dominance 
(associated to Mr = 100 GeV), which leads to larger values for \x — e conversion rates. 





M [GeV] 


(YjY u ) 12 


CRau (current) 


10 i2 


2.7 x 10" 4 


CRxi (current) 


4 x 10 12 


10~ 3 


CR Ti (future: 10~ 16 ) 


2 x 10 10 


5.5 x 10" 6 


CR Ti (future: 10" 18 ) 


2 x 10 9 


5.5 x 10~ 7 



Table 1: Approximate limits on M and (YjY^,)^ from the non-observation of fi — e conversion in 
nuclei. This table includes current experimental bounds and future expected sensitivities [62,63, 

E5J. 

The non-observation of \i — e conversion in gold, yg 7 Au, sets an upper bound on the size of M, 
M ~ 10 12 GeV (the exact value slightly dependent on Mr). 

Conservative limits will be derived from the results obtained for Mr > 1 TeV. We stress that 
lower values for Mr would lead to more stringent limits. This approximate limit on M can be 
translated into constraints for (YjY u )i2, the flavour violating combination of Yukawa couplings 
controlling [i — e transitions. Using Eqs. (|2.6p and ()2.9|) . we find that this limit is 2.7 x 10 -4 . On 
the other hand, /i — e conversion in llTi has a slightly more relaxed experimental bound, which in 

15 This has been studied in [ID] for the non-supersymmetric inverse seesaw case. 
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Figure 5: [i — e conversion rates in 7<j 7 Au (left) and ll^i (right), as a function of M and (YJY u )i2 
for different M R values: M R = 100 GeV (blue), M R = 1 TeV (red) and M R = 10 TeV (black). We 
set Aq = —300 GeV, tan/3 = 10, sign(/x) = + and Bq = 0, and we randomly vary ttlq and Mt/ 2 in 
the range [0,3 TeV]. The horizontal dashed lines represent the current experimental bounds and 
the dotted ones represent the expected future sensitivities. 



turn implies less stringent limits on both M and (YJ Yj,)i2, 4 x 10 12 GeV and 10 -3 , respectively. 
However, there are plans to improve the experimental sensitivities for /1 — e conversion in titanium 
with expected sensitivities in the 10 -18 — 10 -16 range [62,63,65 . Should this be the case, and 
no n — e conversion events be observed, the limits on the parameters of the inverse seesaw model 
would be greatly improved. 

In Table [U we compute approximate limits for M and (YjY u )i2, as obtained from current 
experimental bounds and future expected sensitivities. Such a (rough) estimate is possible due to 
the little dependence of this result on M R and on the cMSSM parameters, as shown in Figure [5j 
However, from these estimates, we can safely conclude that 0(1) Yukawa couplings, as typically 
chosen in the literature, are clearly ruled out since they would induce excessively large cLFV 
effects. This would imply that (YjY u )i2 ~ 1, corresponding to M ~ 10 14 GeV, clearly out of the 
allowed ranged. 



5.4 l t Uj 

Concerning the 3-body decays l{ — > 3£j we focus on \i — > 3e, due to its very challenging experimen- 
tal bound. Since the Z-boson does not differentiate between leptonic flavours, the same behaviour 
and enhancement can be found for the cLFV r 3-body decays into leptons. 

As discussed in [40], the branching ratio for [i — > 3e is greatly enhanced by the Z-penguin. 
In fact, it turns out that this process is only slightly less constraining than fi — e conversion in 
nuclei. In Figure[6l the Br(/x — > 3e) is depicted as a function of M, and for three representative M R 
values, again assuming a degenerate singlet spectrum. As was the case of \i — e conversion in nuclei, 
the dependence on the cMSSM parameters is essentially negligible. Similarly, the scenario with 
M R = 100 GeV is dominated by the non-super symmetric photonic penguins and larger branching 

16 This statement only applies to off-diagonal Yukawa couplings. As explained at the beginning of this section, the 
fix limit, in which the Y v matrices are exactly diagonal, suppresses all cLFV observables and thus 0(V) diagonal 
Yukawa couplings would be allowed. 
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Figure 6: Br(/i — > 3e) function of M and (yJy u ) 12 for different Mr values: M R = 100 GeV 
(blue), Mr = 1 TeV (red) and Mr = 10 TeV (black). The cMSSM parameters are taken as in 
figure The horizontal dashed line represents the current experimental bound. 



ratios can thus be obtained. As can be seen, the corresponding limits on M and (YjY u )i2 would 
be slightly more relaxed than those arising from /i — e conversion in gold and titanium. 

It is straightforward to obtain the following relation between the branching ratios of different 
ti — > 3£j channels 

Bv(£ m ^3£ n ) (YjY u )l n ml m r m ' 

where Tj im are the life-times of the leptons. Equation (j 5 . 3 j) provides a very good approximation to 
the numerical results and allows to conclude that, unless strong cancellations occur in the (YvY v )ij 
terms, the three observables, Br(/z — > 3e), Br(r — > 3e) and Br(T — > 3fi), are expected to lie within 
1 — 2 orders of magnitude. For example, assuming degenerate right-handed neutrinos, a vanishing 
Dirac phase in the neutrino sector, a normal hierarchy and a vanishing lightest neutrino mass 
(m ui = 0), one finds Br(/x — > 3e) ~ 70Br(r — > 3e) ~ 0.7Br(r — > 3/i). An observation of r — > 3e or 
r — > 3fi at rates much larger than those expected for [i — > 3e could only be accommodated with a 
strong cancellation in (YvY v )\2, that would suppress the fj, — e transitions while still allowing for 
r — e and Similar results for this observable in an inverse seesaw framework with an 

extended gauge group have been found in [21]. Finally, in Z-penguin dominated scenarios there 
is a clear correlation between the rates for (i — > 3e and /i — e conversion in nuclei. Numerically, 
we find CR(jU - e,Ti)/Br(/x 3e) ~ 15. 



5.5 r P £i 

We now address r flavour violating decays with a light pseudoscalar meson in the final state, 
r — > P°£i. As done for the previous observables, we present our results for the branching ratios 
as a function of M, and for three different values for Mr. This can be seen in Figure [TJ where we 
focus on the particular case r — > r^i. 

We observe in Figure [7] a clear enhancement which is again due to the dominance of the Z- 
boson mediated contributions in the total amplitude. However the branching ratio remains small 
as the size of the (YJ 1^)23 terms never exceeds O(10 -2 ). 

17 In such a scenario, fi — e conversion in nuclei and fj, — >• e~/ would be suppressed as well. 
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Figure 7: Br(r — > tj/jl) function of M and (y2y v ) 2 3 for different Mr values: Mr = 100 GeV 
(blue), Mr = 1 TeV (red) and Mr = 10 TeV (black). The cMSSM parameters are taken as in 
figure 

In fact, without the Z-penguin contributions, this observable would be much more suppressed 
(by several orders of magnitude). For example, the Higgs- mediated contributions to this process 
were studied in [25| and it was found that 0(1) Yukawa couplings would be required to reach 
observable levels. However, the strong constraints set by [i — e conversion in nuclei preclude this 
possibility, implying that Br(r — > rjfj,) < 10~ 13 . In conclusion, if the inverse seesaw is realised in 
Nature, semileptonic r cLFV decays do not have realistic chances of being observed in the near 
future unless strong cancellations occur in the fj, — e sector. 

Very similar results have been obtained for other mesons, as well as for final states including 
an electron instead of a muon. These observables also exhibit a very little dependence on the 
spectrum: large variations in tuq and Mi/ 2 have a minimal impact on the results for Br(r — > P°£i). 

5.6 ti Iff 

Finally, we discuss our results for 1^ — > ^7. Radiative decays have been intensively studied in 
the literature (see for instance [9,48,68]). However, most of the phenomenological studies have 
focused only on this observable, neglecting the impact of the Z-boson enhancement on the other 
cLFV observables. Therefore, our purpose is to determine how the constraints on the parameters 
of the inverse seesaw extension of the MSSM, derived from the analysis of the Z-boson enhanced 
observables, affect the predictions for the l{ — > rates. In fact, due to the strong limits on 
the Yukawa couplings coming from fi — e conversion in nuclei, Br(£j — > ij'y) is typically below 
the experimental limits, the only exception being the Mr = 100 GeV scenario. This can be seen 
in Figure [8] where, as expected from theoretical arguments, Br(/i — > 3e) > Br(;U —> e^f) in the 
Z-penguin dominated scenario£E 

It is also clear from Figure [8] that for low Mr, the observable Br(/i — > e~f) has very little 
dependence on ijiq and M1/2, whereas for large Mr, one can find very large variations due to the 
different values of the SUSY masses. Similarly, we have checked that the other cMSSM parameters, 
Aq, tan/3, sign(^) and Bq, do not significantly affect the numerical value of Br(/i — > e^f) when 
Mr = 100 GeV, but can become relevant for larger values of Mr. We find that the non-SUSY 

18 This is true for moderate values of mo and M-1/2 or , in other words, when the superparticles running in the 
loop are reasonably heavy. Similarly, low Mr spoils this feature by increasing the non-SUSY contributions. 
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Figure 8: Br(/x — > e7) function of M and {YJY U ) 12 for different Mr values: Mr = 100 GeV 
(blue), Mr = 1 TeV (red) and Mr = 10 TeV (black). The cMSSM parameters are taken as in 
figure The horizontal dashed line represents the current experimental bound. 



contributions to Br(/i — > ej) become relevant only for Mr < 1 TeV, and in fact, for Mr = 100 GeV 
the non-SUSY contributions totally dominate, so that all dependence on mo, Mij 2 and on the rest 
of cMSSM parameters disappears. These results are in agreement with those found in [9], where 
it was shown that non-SUSY contributions can enhance Br(£j — > £jj) in the inverse seesaw if the 
singlet fermions are light and have relatively large mixings with the active neutrinos. This result 
has been confirmed by our numerical study, where we found that Br(^j — > ij^y) could be enhanced 
by some orders of magnitude when Mr <C m\y. One can thus find regions of the Mr — M plane 
where Br(^ — > 3e) < Br^ — > e'j), contrary to what occurs when SUSY contributions dominate in 
both observablea 19 !. 



5.7 Other observables and benchmark points 

To conclude the numerical analysis, we address some aspects not fully covered in the previous 
sections. We also summarise the most relevant features and results using some specific benchmark 
points. 

Concerning flavour violating neutral boson decays, some additional comments are in order. 
The cLFV decays of the Z-boson, Z — > with i ^ j, are not enhanced^! and the corresponding 
branching ratios lie below experimental reach [87] . Similarly, the branching ratios for H° — > £i£j 
are very small. This process is known to have a non-decoupling behaviour leading to non-neglible 
contributions even for very heavy sparticles [44] . However, in the present case, the small Yukawa 
couplings strongly suppress the corresponding rate. We have also studied collider observables such 
as Br(x° — > Xi^i&j) an d slepton mass splittings, Am,-. We have found that they have very little 
relevance in this model as expected from the previous discussion in Section T4.6I 

19 The Mr — M plane is not only constrained by cLFV observables, but also by Non-Standard Interactions 
(NSI) |85J . In this model, for degenerate singlets, these bounds can be translated into fix > 20 eV. This in turn 
implies that, for a given value of M, one can always find a lower bound on Mr. For example, for M = 10 12 GeV 
the NSI bound implies Mr > 135 GeV and thus the blue dots in Figure[H]with M > 10 12 GeV, which are ruled out 
by jj, — > e7, are also disfavoured by NSI. For a non-degenerate spectrum the previous estimate does not hold. 

20 This process has been studied in detail in [86], where the relevant expressions are given. Recently, good 
agreement, up to an overall factor of 2 in the wave function contributions, was found in |41j . 
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Figure 9: R^ e -y (blue) and Rau (red) as a function of M. The Mr matrix is given by Mr = 
diag(Afja, lOTeV, lOTeV), with M R1 = 30 GeV (solid lines), M B1 = 100 GeV (dashed lines) and 
Mri = 1 TeV (dotted lines). The cMSSM parameters take the values m = 500 GeV, M ly / 2 = 1 
TeV, A = -300 GeV, B = and tan /3 = 10. 



For completeness, we briefly discuss the case of non-degenerate right-handed neutrinos, Mr = 
diag(M_Ri, Mr2, Mrs), with Mri ^ Mrj. In our analysis we considered three different values for 
Mm (30 GeV, 100 GeV and 1 TeV), fixed M R2 = M R3 = 10 TeV and varied M in the range 
[10 8 , 10 13 GeV]. In order to allow for a comparison with the degenerate case (setting, for example, 
Mri = Mr2 = Mrs = 10 TeV) and to identify the effect of the right-handed neutrino spectrum 
on the cLFV observables, we have taken a common M for the three right-handed neutrinos by 
adjusting the corresponding fixi- This allows to have identical values for the Yukawa couplings in 
the degenerate and non-degenerate cases. To better evaluate the difference between the degenerate 
and non-degenerate cases, we further define the ratios 



BrQz -> e7) non _ deg _ CR(/i - e, Au) non _ deg 

Br(/^e 7 ) deg ' * Au " CR(/i-e,Au) deg 



The results can be seen in Figure [9] where the cMSSM parameters take the arbitrary (but repre- 
sentative) values m = 500 GeV, M 1/2 = 1 TeV, A = -300 GeV, B = and tan/? = 10. As 
expected, Br(/i — > ej) is affected by a change in the singlet spectrum. However, one can see that 
the relative difference (for this particular choice of parameters) can reach at most ~ 10%. The 
change in the Z-boson mediated processes is negligible: for example, the relative change in the 
fi — e conversion rate in gold is always below 2%, independently of the value of M. This is another 
example of the little impact that the mass spectrum has on the Z-boson mediated processes. In 
conclusion, a non-degenerate right-handed neutrino spectrum can indeed induce changes to the 
observables, but these are typically small and the essence of the results derived here is unaffected 
by the nature of the right-handed neutrino spectrum^!. 

Finally, we present our results for the cLFV observables for the two benchmark points of 
Table [2j Point A is associated to a Higgs mass m^o = 126.5 GeV, in agreement with the latest 
LHC [761 [77] and Tevatron [78] results© We note that, although the inverse seesaw contains 
additional degrees of freedom beyond those of the cMSSM, the smallness of the Yukawa couplings 
implies that any correction to the Higgs mass, m h o, arising from the singlet sector will be very 



21 There is an exception to this general statement: in the non-degenerate case, the R matrix is indeed relevant 
even if it is real [46]. 

22 Our numerical routines for the computation of the Higgs mass, based on SPheno 3.1.9 69,70], agree with the 
results of [79H81] . within the usual theoretical uncertainty of ~ 1 — 2 GeV. 
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Table 2: Benchmark points used in the numerical analysis. Point A leads to a Higgs mass of 
m(/i°) = 126.5 GeV, in accordance with the latest LHC results. Point B is an example of a very 
heavy supersymmetric spectrum. 

smal@. Point B exemplifies a very heavy supersymmetric spectrum with associated gluino and 
squark masses larger than 6 TeV and lightest chargino and sneutrino masses of 2.3 TeV and 3 
TeV, respectively. It has been selected in order to show explicitly the non-decoupling behaviour 
of the Z-mediated processes. 

The results for the cLFV observables are given in Tableland, for completeness, we display in 
Table H] the current experimental bounds and future sensitivities for the same cLFV observables. 
For each point we have taken three different spectra in the singlet sector: (1) M R \ = Mr2 = 
M m = 100 GeV, (2) M R1 = M R2 = M R3 = 1 TeV, and (3) M R1 = 50 GeV, M R2 = 500 
GeV, M R 3 = 1 TeV. The first case serves to show the influence of the non-SUSY contributions 
to £i — > Ij^f, the second one is a standard degenerate scenario where the SUSY contributions 
dominate all processes, while the third case is a non-degenerate scenario with an intermediate 
situation. In all cases we fix M = 4 x 10 11 GeV. 

Firstly, we observe that Br(/j — > £jj) shows very little change when going from scenario A- 
1 (with relatively light sparticles) to B-l (with very heavy sparticles). As discussed in Section 
15.61 in scenarios with very low masses for the right-handed neutrinos, the main contributions to 
Br(/j — > £j-f) are non-supersymmetric. Therefore, they are obviously unaffected by the size of the 
SUSY parameters. We also point out that, although the right-handed neutrinos are light in these 
two scenarios, the bounds from NSI [85] are satisfied due to the smallness of the Yukawa couplings. 
Another interesting remark that can be made in association to Table [3] is the different behaviour 
of the radiative £{ — > £jj and the Z- mediated processes when going from the A-2 point to the 
B-2 point. While Br(£j — > £j-f) is clearly reduced, in agreement with the well-known dependence 
on the SUSY spectrum, oc m-g^gy, the rates for the Z- mediated processes hardly change. This 
is a clear indication of the non-decoupling behaviour due to the Z-penguins. Finally, points A-3 
and B-3 show an intermediate case, with one light right-handed neutrino, yielding results quite 
similar to those in points A-2 and B-2. Thus, as discussed before, the nature of the right-handed 
neutrinos spectrum is irrelevant. 

6 Conclusion 

The supersymmetric inverse seesaw is a very attractive extension of the MSSM, with neutrino 
masses generated by TeV-scale mediators, in association with potentially large Yukawa couplings. 
In this work, we have studied in detail the predictions for several lepton flavour violating observ- 
ables, focusing on those mediated by Z- boson exchange. These are particularly interesting when 
the non-supersymmetric contributions are small. 

We have found that, due to the non-decoupling behaviour, the Z-penguins totally dominate 
the cLFV amplitudes in most of the parameter space, especially in scenarios where the right- 
handed neutrinos and the supersymmetric particles have masses larger than ~ 500 GeV. In those 

23 For large Yukawa coupling, an enhancement to the Higgs mass can be found. However, this is often in conflict 
with bounds from LFV l88l . 
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Table 3: Results for several lepton flavour violating processes for the benchmark points of Table [2] 
supplemented with three different singlet spectra: (1) Mri = Mr2 = Mr3 = 100 GeV, (2) 
M R1 = M R2 = M m = 1 TeV, and (3) M R1 = 50 GeV, M R2 = 500 GeV, M m = 1 TeV. 
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Table 4: The current experimental bounds and future sensitivities for the cLFV observables. 

cases, [i — e conversion in nuclei is the most constraining observable, clearly more restrictive than 
fx —> ej. As a result of the large enhancement provided by the Z-penguins, one can set very strong 
constraints on the size of the flavour violating couplings. 

Finally, we emphasize that the bounds obtained in this work apply to the off-diagonal elements 
of the Yukawa couplings. One can also set bounds on the diagonal ones using the invisible width 
of the Z boson. This effect is negligible in our case due to the smallness of the Yukawa couplings. 
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